Observation of an Anomalous Number of Dimuon Events in a High Energy Neutrino 
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A search for long-lived neutral particles (7V°'s) with masses above 2.2 GeV/c 2 that decay into at 
least one muon has been performed using an instrumented decay channel at the NuTeV experiment 
at Fermilab. Data were examined for particles decaying into the final states fifi, fie, and [im. Three 
jx\x events were observed over an expected Standard Model background of 0.069 ± 0.010 events; no 
events were observed in the other modes. 

. PACS numbers: 14.80.-j, 13.85.ftm, 12.60. Jv, 13.15+g 



A number of extensions to the Standard Model (SM) 
of elementary particle physics predict new particles with 
reduced couplings to normal quarks and leptons. Some 
of these particles, which we shall refer to as N°'s, may 
have zero electric charge, long lifetimes, and small in- 
teraction rates with ordinary matter. Examples include 
neutral heavy leptons (NHLs) or heavy sterile neutri- 
nos and neutral supersymmetric particles Q. The 
N° may be either pair-produced in hadronic interactions 
or occur as a daughter particle in weak decays of mesons 
through mixing with SM neutrinos. Decays of the N° can 
proceed through M /± or Z° mediated decays via mixing, 
or through iZ-parity violating supersymmetric processes. 
We report here the results of a search for N° particles 
in the mass region above 2.2 GeV/c 2 that decay into fi- 
nal states with at least one muon and one other charged 
particle, using the NuTeV detector at Fermilab (E815). 
This analysis extends our results for a previous search for 
JV°'s in the mass region 0.3 to 3.0 GeV/c 2 with at least 
one final state muon ||, and for a search for 7V°'s with 
mass < 0.3 GeV/c 2 that decay to e + e~^ ||. 

NuTeV received 2.54 x 10 18 800 GeV/c protons during 
the 1997 fixed-target run at Fermilab with the detector 
configured for this search. The proton beam was inci- 
dent on a one-interaction-length beryllium oxide target 
at a targeting angle of 7.8 mr with respect to the detec- 
tor. A sign-selected quadrupole train (SSQT) |t| focused 



either positive (45% of the running) or negative (55% of 
the running) secondary n and K mesons into a 440 m 
evacuated pipe pointed toward the NuTeV experiment. 
Surviving neutrinos, and possibly JV°'s, traversed ^850 
meters of earth-steel shielding before reaching the NuTeV 
decay channel. 

The decay channel region (Fig. |Tj), located 1.4 
km downstream of the production target, was de- 
signed to search for exotic neutral particle decays in a 
(1.27 x 1.27 x 34) m 3 fiducial volume. The volume con- 
sisted of large helium-filled plastic bags interspersed with 
drift chambers (DC). An array of plastic scintillation 
counters at the upstream end of the decay channel pro- 
vided a means to veto charged particles produced in the 
upstream shielding. The Lab E neutrino detector 
downstream of the decay channel provided final state 
particle energy measurement and identification. This de- 
tector consisted of a 690 ton iron-scintillator sampling 
calorimeter, instrumented with DC for charged particle 
tracking, followed by a toroid muon spectrometer. Sets 
of hits in the calorimeter DC were linked to tracks found 
with the decay channel DC system to determine parti- 
cle identification. Muons were identified by penetration, 
electrons by their characteristic short clusters of hits, 
and charged hadrons by their elongated hit patterns. All 
charged hadrons were assumed to be pions. The particle 
identification algorithm was tuned with electrons, pions, 
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and muons provided by the NuTeV calibration beam |g| . 
The probability to mis-reconstruct a dimuon event as a 
fie or fiit event was 0.2%; /xe and [iir events were mis- 
reconstructed as \x\x events with a probability of 0.4% 
and 1.3% respectively. Muon energy was determined by 
either spectrometer measurement {o v jp = 11%), range in 
the calorimeter (er p = 310 MeV), or multiple scattering 
(<T p /p — 42% at p = 50 GeV/c). Only muons penetrat- 
ing to the toroid spectrometer could be charge-identified. 
The hadronic energy resolution of the calorimeter was 
a/E= (0.024 ±0.001) © (0.874 ± 0.003)/^E(GeV); the 
electromagnetic energy resolution, cr/E = (0.04±0.001)© 
(0.52 ± 0.01)/V£(GeV) §. 

Event selection criteria were developed to minimize 
known backgrounds while maintaining efficiency for a 
possible N° signal. A series of cuts isolated events with 
exactly two well-reconstructed tracks forming a vertex 
within the decay channel fiducial volume and having no 
charged particle identified in the upstream veto system. 
Both tracks were required to be well-reconstructed and 
have an associated calorimeter cluster, with at least one 
of the tracks identified as a muon. The track and vertex 
quality criteria were numerically the same as in Ref . || . 
The probability for a signal event to fail this is <2%. 
The vertex position was required to be within the de- 
tector fiducial volume; in addition, the longitudinal dis- 
tance from the vertex position to any drift chamber was 
required to be greater than the larger of 101.6 cm and 
3er z , with a z the longitudinal vertex position error. A 
third track which formed a downstream vertex with one 
of the two initial tracks was permitted, to allow for S- 
ray emission. Cosmic ray tracks (which are suppressed 
by the fast gate timing of the neutrino beam) were re- 
moved by requiring the slope of each track relative to the 
beam direction be less than 100 mr. Muons, hadrons, 
and electrons were required to have an energy greater 
than 2.2 GeV, 10 GeV, and 10 GeV, respectively, with 
an additional total energy cut of 12 GeV applied to fifi 
events. In order to isolate high mass events, a trans- 
verse mass cut mj- > 2.2 GeV/c 2 was applied, with 
niT = \Pr\ + \/Pt~+ m 'v ' Pt the component of the to- 
tal reconstructed momentum perpendicular to the beam 
direction, and my the invariant mass of the visible par- 
ticles. 

Further "clean event" cuts were applied to reduce the 
level of the dominant deep-inelastic neutrino scattering 
(DIS) backgrounds. DIS events typically had large track 
multiplicities, many drift chamber hits, and extra unas- 
sociated clusters in the calorimeter. Clean cuts required: 
(1) three or fewer tracks in any one DC view, (2) three 
or fewer DC hits in any view of the first chamber down- 
stream of the vertex, (3) at least one DC view with 
fewer than eight hits total in the first two chambers 



downstream of the vertex, (4) no energy clusters in the 
calorimeter not associated with tracks, and (5) no tracks 
identified as electrons with missing hits in either view of 
the first two chambers downstream of the vertex. The fi- 
nal cut removed events where a photon from the primary 
vertex converted to e + e~ and was reconstructed as an 
electron. 

Detailed MC simulations of physics processes includ- 
ing detector effects were used to estimate possible back- 
grounds to the N° signal. Two major classes of physics 
processes considered included DIS pOj l, r esonance pro- 
duction [fi"l| , and diffractive scattering |l2[] by neutrinos; 
and decays and interactions of hadrons and photons pro- 
duced in neutrino interactions. Particular attention was 
given to known sources of dimuon or dimuon-like pro- 
duction: DIS, resonance, and diffractive production of 
charm; neutrino trident production; fi + fi~ vector meson 
decays; electromagnetic muon pair production; low mul- 
tiplicity DIS accompanied by a secondary pion or kaon 
decay; and decays of mesons produced by neutrino 
interactions in the decay channel or surrounding mate- 
rial. Approximately 10 3 times as many neutrino interac- 
tions as in the data were simulated in the decay channel 
volume; in addition, a large sample of DIS events was 
generated in the material surrounding the decay channel. 
Event generators fed a GEANT-based Jl3| detector sim- 
ulation that produced hit-level simulations of raw data 
including DC inefficiencies and noise. MC events were 
processed using the same analysis routines used for the 
data. A number of other possible sources were studied 
as well and found to be negligible. 

Background calculations were normalized to data us- 
ing charged-current DIS interactions in the decay channel 
DC. Events in this sample were required to pass the fol- 
lowing five normalization cuts: a vertex within the trans- 
verse fiducial volume ; a z vertex within 76.2 cm of a DC; 
no upstream veto; >1 GeV energy deposit in the front of 
the calorimeter; and one toroid-analyzed muon matched 
to a decay channel track. The MC was normalized to 
match the total number of data events with two or more 
tracks; the error on this normalization is 9%. Figure || 
shows a comparison of the longitudinal vertex resolution 
between data and MC with the z vertex requirement re- 
moved. 

MC events were also compared to another data control 
sample as a check on the quality of the simulation. For 
this sample, the vertex was required to be within the de- 
cay channel transverse fiducial volume but the z position 
was allowed to be either in the chambers or the helium. 
Tight track angle cuts were imposed, and there was a 
strict requirement on veto system activity. The majority 
of these events were from interactions in the chamber ma- 
terial or from interactions in the laboratory floor. Of 495 
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events in the data, 116 had vertices reconstructed in the 
helium at least 101.6 cm from the nearest DC. This can 
be compared to the MC, which predicted (514 ±82) total 
events and (96 ± 15) events reconstructed in the helium. 
Because loose vertex quality requirements allowed mis- 
reconstructed interactions in the chambers and floor to 
enter this sample, only 35% of the vertices reconstructed 
in the helium were actually due to v-He interactions. 

After all cuts, the expected background to the N° 
search was 0.069 ± 0.010 events in fifj, mode, 0.13 ± 0.02 
events in fj,e mode and 0.14 ± 0.02 events in /itt mode. 
Table | summarizes the background components for /i/i 
mode; DIS clearly dominates. 

Before looking at the data in the signal region, we per- 
formed a series of analyses on other fiducial and kine- 
matic ranges. These included using: (1) identical anal- 
ysis cuts applied to events within 15.2 cm of a DC (the 
"chamber region"); (2) the chamber region with loosened 
cuts to increase fin acceptance; (3) the "intermediate re- 
gion" between 15.2 and 101.6 cm from the chambers, with 
otherwise standard analysis cuts; and (4) events with 
well-reconstructed two-track vertices where the tracks 
were both identified as pions. Measurements agreed with 
MC predictions within 1.5cr in all cases. For example in 
the "chamber region" sample (1), two two-track events 
were observed in the data with 3.5 predicted by the MC. 

Upon examining the signal region, three mi(y) events 
were observed, which is considerably above the predicted 
background. No /ie or events were observed, which 
is consistent with background estimates. Table || sum- 
marizes event reconstruction characteristics of the 
events. The vertex and track probabilities have been 
measured with both the signal and background Monte 
Carlos and represent the probability that the vertex (or 
track) would reconstruct with a similar (or worse) qual- 
ity- 
Observation of three fifi data events prompted further 
tests comparing data to MC predictions with relaxed 
cuts. These included gradually removing the cuts and 
alternately releasing and then returning individual cuts 
and cut-pairs. These studies revealed no significant dis- 
crepancy between data and the background calculations 
(except for the original 3 fifi events). 

The three events have some features consistent with 
a N° decay hypothesis. The events pass the analysis cuts, 
where the background is estimated to be 0.069 events. 
All three occur well within the fiducial volume away from 
the chambers and are evenly distributed throughout the 
decay channel. The transverse mass, invariant mass, and 
missing Pt are all consistent with the decay of an TV 
with a mass of about 5 GeV/c 2 . Since only 1(0) fie and 
2(0) fiTr events (consistent with MC expectations) were 
observed in the chamber (helium) data, it is unlikely that 



Mi events are related to low multiplicity neutrino events 
followed by 7r and K decay. 

Unlike the background, in both an NHL or neutralino 
model one would expect the /Z7r rate to be highly sup- 
pressed relative to leptonic decays. However, for a 5 
GeV/c 2 NHL model, one would expect 1.4 times more 
\iev events ||. A neutralino model, on the other hand, 
can accommodate the observation of either only [i\x or a 
combination of fifj, and fie candidates by selecting appro- 
priate couplings. 

However, several aspects of the candidate events are 
similar to those from neutrino interaction backgrounds, 
and might be indicative of unaccounted-for sources or 
a statistical fluctuation. Globally, the events share one 
feature that is improbable for an iV° decay hypothe- 
sis. All three events have a muon energy asymmetry 
A > 0.83, where A = (\E t - E 2 \) / (E 1 + E 2 ). For DIS 
background, the probability for three \x\i events which 
pass the signal cuts to have the observed energy asym- 
metry is 38%. The probability that this occurs in a weak 
decay hypothesis (lj] is less than 0.5% (including accep- 
tance). All three events occurred during the higher rate 
i'- mode rather than F-mode running periods. In the two 
events where the charge of the higher-momentum muon 
can be measured, it has the same sign as expected for the 
leading muon in a charged-current neutrino interaction. 
Event kinematics (Mj, M^, Pt) are also consistent with 
DIS characteristics. However, the observed number of 
events is inconsistent with expected neutrino interaction 
background. 

In summary, NuTeV has observed three \i\x events, zero 
fxiT, and zero /ie events with transverse mass above 2.2 
GeV/c 2 . The expected backgrounds were 0.069 ± 0.010, 
0.14 ±0.02, and 0.13 ±0.02 events, respectively. The rate 
corresponding to the observed three events is not consis- 
tent with Standard Model processes we have identified 
and the source of the events is not clear. NuTeV is in- 
sensitive to NHL production in this mass region, but can 
set an interesting limit on neutralino production (deter- 
mined by calculating one-sided limits using a frequentist 
approach without background subtraction |^5|). NuTeV 
is the first experiment to set limits on the production 
of long-lived neutralinos in this mass range which de- 
cay by i?-parity violation. This limit (Fig. |^), although 
motivated by a neutralino hypothesis, is a generic limit 
applicable for any model of neutral particle production 
at the target Q. 
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FIG. 1. Schematic of the NuTeV decay channel with ex- 
ample N° — > /iir decay. 
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FIG. 2. Distance of longitudinal vertex position from the 
closest chamber for events in the normalization sample with- 
out the z vertex requirement. (Crosses: data; histogram: 
MC; arrows: values for the three observed mjl events) 
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FIG. 3. NuTeV limit on neutralino production. The plot 
gives the limit on the differential p-p cross section to produce 
an N° versus the decay length of the N° in the laboratory 
coordinate system. This limit is generic for an TV produced 
at the target. 
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TABLE I. Estimated rates of background to the 
TV — > mi{v) search 



Source MM(^) events 

DIS events (6.8 ± 1.0) x 1(T 2 

Diffractivc charm (1.3 ± 0.1) x 10~ 3 

Diffractive tt (1.9 ± 0.1) x 10~ 4 

Diffractive K (4.0 ± 0.3) x 10~ 7 

K° L decays from shielding (3.9 ± 3.9) x 10~ 4 

Other sources < 2.5 x 10~ 4 

Total fifj,(i/) Background (6.9 ± 1.0) x 10~ 2 



TABLE II. Kinematic and reconstruction quantities for the 
three candidate N° — > fJ,fJ,(v) events. The variables refer to 
the muons energies (E M i, E M 2), missing transverse momentum 
(-Prmiss), two muon invariant mass (m; nv ), transverse mass 
(niT, see text), transverse vertex position (v x ,v y ), longitudi- 
nal distance to nearest drift chamber (|Az|), vertex proba- 
bility (Vvert), first muon track probability ("P M i), and second 
muon track probability (V^)- The probabilities are measured 
the with signal(background) Monte Carlos. The sign on the 
muon energy refers to the charge of the muon (if measured) . 



Event 1 2 3 



E M i (GeV) 


-77.4 


-92.0 


±48.0 


E M2 (GeV) 


±2.56 


±5.85 


±4.34 


P T m iss (GeV/c) 


2.42 


1.41 


2.07 


m inv (GeV/c 2 ) 


1.10 


0.88 


3.57 


m T (GeV/c 2 ) 


5.08 


3.08 


4.66 


v x (cm) 


-46.5 


48.0 


-57.5 


v y (cm) 


3.4 


-38.3 


15.0 


\Az\ (cm) 


193 


657 


-188 


invert) 


0.81(0.94) 


0.043(0.48) 


0.011(0.30) 


(V,i) 


0.95(0.96) 


0.034(0.18) 


0.51(0.73) 


(V»2) 


0.49(0.70) 


0.034(0.18) 


0.15(0.39) 



